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Abstract. We discuss a phenomenological approach to the description of unstable 
f"|^ ' vehicle motion on multilane highways that could explain in a simple way such 

O ' observed self-organizing phenomena as the sequence of the phase transitions "free 

^ , flow —^ synchronized motion —^ jam" and the hysteresis in them. 

We introduce a new variable called order parameter that accounts for possible 

correlations in the vehicle motion at different lanes. So, it is principally due to 

"many-body" effects in the car interaction in contrast to such variables as the 

— L ' mean car density and velocity being actually the zeroth and first moments of the 

(-H , "one-particle" distribution function. Therefore, we regard the order parameter as an 

Q ■ additional independent state variable of traffic flow and formulate the corresponding 

Q_) ' evolution equation governing the lane changing rate. 

In this context we analyze the instability of homogeneous traffic flow manifesting 
itself in both of these phase transitions and endowing them with the hysteresis. 
wC [ Besides, the jam state is characterized by the vehicle flows at different lanes being 

VQ . independent of one another. 



1 Introduction 



The existence of a new basic phase in vehicle flow on multilane highways 
called the synchronized motion was recently discovered by Kerner and Re- 
j^ ' hborn ll|, impacting significantly the physics of traffics as a whole. In particu- 

lar, it turns out that the spontaneous formation of moving jams on highways 
proceeds mainly through a sequence of two transitions: "free flow — > syn- 

■ O I chronized motion -^ stop-and-go pattern" |p|. Besides, all these transitions 

exhibit the hysteresis P0,H. As follows from the experimental data 0,3, i| 
fj I the synchronized mode is essentially a multilane effect. Recently Kerner |5,3| 

assumed that the transition "free flow -^ synchronized mode" is caused by 
"Z"-like form of the overtaking probability depending on the car density. 

j^ ■ There have been proposed several macroscopic models dealing with mul- 

tilane traffic flow 0J|,|,^jll||l|,|l3|,0 . Both these models specify the traffic 



dynamics completely in terms of the car density p, mean velocity v, and, 
may be, the velocity variance 9 or ascribe these quantities to the vehicle 
flow at each lane individually. Nevertheless, a quantitative description of the 
synchronized mode is far from being developed well because of its complex 
structure |5|,|6|]. In particular, it can form the totally homogeneous (i) and 



homogeneous-in-speed (ii) flows m. Especially in the latter case there is no 
explicit relationship between the mean car velocity v and density p, with the 
value of V being actually constant and less then that of free flow. The other 
important feature is the key role of some cars bunched together and traveling 
much faster than the typical ones, which enables to regard them as a special 
car group ||l|. Therefore, in the synchronized mode the function of car distri- 
bution in the velocity space should have two maxima and we will call such 
fast car groups platoons in speed. These features of the synchronized mode 
have been substantiated also in |l5[ using single-car-data. In particular, it has 
been demonstrated that the synchronized mode exhibits small correlations 
between fluctuations in the car flow, velocity and density. There is only a 
strong correlation between the velocities at different lanes taken at the same 
time and decreasing sufficiently fast as the time difference increases. By con- 
trast, there are strong long-time correlations between the flow and density in 
the free flow state as well as the stop-and-go mode. 

Keeping in mind a certain analogy with aggregation processes in physical 
systems Mahnke et al. ]l^,|l3 proposed a kinetic model for the formation of 
the synchronized mode treated as the motion of a large car cluster. In the 
present paper following practically the spirit of the Landau theory of phase 
transitions we develop a phenomenological approach to the description of 
this process. We ascribe to the vehicle flow an additional internal parameter 
will be called below the order parameter h e (0, 1) characterizing the possible 
correlations in the vehicle motion at different lanes and write for it a governing 
equation. For the car motion where drivers do not change lane at all we set 
h — 0, in the opposite limit h = 1. 

2 Order Parameter and the Individual Driver Behavior 

For fixed values of p and v the order parameter h is assumed to be uniquely 
determined, thus, for a uniform vehicle fiow we write: 

dh ,, ^ 

r— = -<P{h,p,v) , (1) 

where r is the delay time and the function <P{h, p, v) fulfills the inequality: 

We note that the time r characterizes the delay in the driver decision of 
changing lanes but not in the control over the headway, so, this delay can be 
prolonged. The particular value h(v,p) of the order parameter results from 
the compromise between the danger of an accident during changing lanes 
and the will of driver to move as fast as possible. Obviously, the lower is the 
mean vehicle velocity v for a fixed value of p, the weaker is the lane-changing 
danger and the stronger is the will to move faster. Besides, the higher is the 



vehicle density p for a fixed value of w, the stronger is this danger (here the 
will has no effect). Thus, the dependence h{v,p) is an decreasing function of 
V and p, so, due to (||): 

^>'' ^p>'' (^) 

with the latter inequality being caused by the danger effect only. Equation 
(El) describes actually the behavior of the drivers that prefer to move faster 
than the statistically mean vehicle and whose readiness for risk is greatest. 
Exactly this group of drivers (platoons in speed) govern the value of h. 

There is, however, another characteristics of the driver behavior, it is the 
mean velocity v = '&{h, p) chosen by the statistically mean driver taking into 
account also the danger resulting from the frequent lane changes by the "fast" 
drivers. Following typical assumptions the velocity ^{h, p) as a function of p 
is considered to be decreasing: 

d-d 

— — < and p^{p) -^ a.s p ^ po , (4) 

where po is the upper limit vehicle density on road. In general, the dependence 
of t9(/i, p) on h should be increasing for small values of the vehicle density, 
/9 <C poi because in this case the lane-changing makes no substantial danger 
to traffic and practically all the drives can pass by vehicles moving at lower 
speed without risk. By contrast, when the vehicle density is sufficiently high, 
p '^ Po, the lane-changing is due to the car motion of the most "impatient" 
drivers whose behavior makes an additional danger to the main part of other 
drivers and the velocity i?(/i, p) has to decrease as the order parameter h 
increases. For certain intermediate values of the vehicle density, p ~ pc, this 
dependence is to be weak as well as near the boundary points, so: 

dd dd d-d 

— > for p < pc , ^ < for p > pc , ^ = at /i = 0, 1 . (5) 

Then the governing equation (nl) takes the form: 

T— = -(j){h, p) , where (j){h, p) =^ ^[h, p, d{h, p)] (6) 

at 

and the condition 0(/i, p) = specifies the steady state dependence h{p) of 
the order parameter on the vehicle density. 

Let us, now, study properties and stability of this steady state solution. 
From Eq. (^) we get 

d(j) d<P d-P dd d(j) d^ d<P dd 

—— — 1 , — - = \ . (7) 

dh dh dv dh ' dp dp dv dp 

As mentioned above, the value of d<P/dp is solely due to the danger during 
changing lanes, so this term can be ignored until the vehicle density p be- 
comes sufficiently high. Thus, in a certain region p < ph < po the derivative 
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Fig. 1. The region of the traffic flow instabihty in the hp-plane and the form 
of the curve h{p) displaying the dependence of the order parameter on the 
vehicle density. 



d(j)/dp ~ {d<P / dv){d'd / d p) < by virtue of (||) and (Q) and the function h{p) 
is increasing or decreasing for d(j)/dh > or d(t)/dh < 0, respectively. This 
statement follows directly from the relation dh/dp = — {d(j)/dp) {d(f)/dh)~ . 
For long-wave perturbations of the vehicle distribution on a highway the 
density p can be treated as a constant. So, according to the governing equa- 
tion (0), the steady-state traffic flow is unstable if dcfy/dh < 0. Due to (0) 
and (H) the flrst term in the expression for d(j)/dh in (M) is dominant in the 
vicinity of the lines /i = and h = 1, thus, in these regions the curve h{p) 
is increasing and the steady state traffic flow is stable. For p < pc the value 
d'd/dh > 0, inequality (||), and, thereby, the region {0 < h < 1, < p < pc} 
corresponds to the stable vehicle motion. However, for p > pc there can be 
an interval of the order parameter h where the derivative d(j)/dh changes the 
sign and the vehicle motion becomes unstable. Therefore, as the car density 
p grows causing the increase of the order parameter h it can go into the in- 
stability region wherein dh/dp < 0. Under these conditions the curve h{p) is 
to look like "S" (Fig. |^a) and its decreasing branch corresponds to the un- 
stable vehicle flow. The lower increasing branch matches the free-flow state, 
whereas the upper one should be related to the synchronized phase because 
it is characterized by the order parameter coming to unity. 



3 Phase Transitions and the Fundamental Diagram 

The obtained dependence h{p) actually describes the first order phase tran- 
sition in the vehicle motion. Indeed, when increasing the car density exceeds 
the value pi the free flow becomes absolutely unstable and the synchronized 
mode forms through a sharp jump in the order parameter. If, however, after 
that the car density decreases the synchronized mode will persist until the 
car density attains the value P2 < Pi- It is a typical hysteresis and the re- 
gion {p2,Pi) corresponds to the metastable phases of traffic flow. It should 
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Fig. 2. The mean vehicle velocity (a) and the vehicle flux (b) vs. the vehicle 
density for the limit values of the order parameter /i = 0, 1 as well as the 
resulting fundamental diagram (c). 



be noted that the stated approach to the description of the phase transi- 
tion "free flow -^ synchronized mode" is rather similar to the hypothesis by 
Kerner 0,0 about "Z"-like dependence of the overtaking probability on the 
car density that can cause this phase transition. 

Let us, now, discuss a possible form of the fundamental diagram showing 
2 = p'd\(^ where, by definition, i?[p] = ^[/i(p),p]. Fig. ||a displays the de- 
pendence 'd{h, p) of the mean vehicle velocity on the density p for the fixed 
limit values of the order parameter ft, = or 1. For small values of p these 
curves practically coincide with each other. As the vehicle density p grows 
and until it comes close to the critical value pc where the lane change danger 
becomes substantial, the velocity i9(l, p) practically does not depend on p. So 
at the point pc at which the curves ^{l,p) and d{0,p) meet each other the 
former curve, 'd{l,p), is to exhibit sufficiently sharp decrease in comparison 
with the latter one. Therefore, on one hand, the function ji{p) = pi?(l, p) has 
to be decreasing for p > pc- On the other hand, at the point pc for h -^ \ 
the effect of the lane change danger is not extremely strong, it only makes 
the lane change ineffective, d-d/dh « (compare (g|)). So it is reasonable 
to assume the function jo(p) = P'&i^TP) increasing neat the point pc- Under 
the adopted assumptions the relative arrangement of the curves jo (p) , ji (p) 
is demonstrated in Fig. 0b, and Fig. 0c shows the fundamental diagram of 
traffic flow resulting from Fig. ^ and Fig. 0b. 

The developed model predicts also the same type phase transition for 
large values of the order parameter. In fact, in an extremely dense traffic flow 
changing lanes is sufficiently dangerous and the function ^(ft, w, p) describing 
the driver behavior is to depend strongly on the vehicle density as p — > po- 
In addition, the vehicle motion becomes slow. Under such conditions the 
former term in the expression for d^jdp in (jj) should be dominant and, so, 
d(j)/dp > and the stable vehicle motion corresponding to d(t>/dh > is 
characterized by the decreasing dependence of the order parameter h{p) on 



the vehicle density p for p > ph- Therefore, as the vehicle density p increases 
the curve h{p) can again go into the instability region (in the Zip-plane), 
which has to give rise to a jump from the synchronized mode to a jam. 
The latter matches small values of the order parameter h (Fig. Hb), so, it 
should comprise the vehicle flows along different lane where lane changing is 
depressed, making them practically independent of one another. 

4 Conclusion 

We have introduced an additional state variable of the traffic flow, the order 
parameter h, that accounts for internal correlations in the vehicle motion 
caused by the lane changing. Since such correlations are due to the "many- 
body" effects in the car interaction the order parameter is regarded as an 
independent state variable. Keeping in mind general properties of the driver 
behavior we have written the governing equation for this variable. 

It turns out that in this way such characteristic properties of the traffic 
flow instability as the sequence of the phase transitions "free flow — » syn- 
chronized motion — > jam" can be described without additional assumptions. 
Moreover, in this model both the phase transitions are of the first order and 
exhibits hysteresis. Besides, the synchronized mode corresponds to highly cor- 
related vehicle flows along different lanes, /i w 1, whereas in the free flow and 
the jam these correlations are depressed, h <^1. So, the jam phase actually 
comprises mutually independent car flows along different lanes. 
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